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We present a simple microscopic model to show how fluctuating two-level systems in a Josephson 
junction tunnel barrier of thickness L can modify the potential energy of the barrier and produce 
critical current noise spectra. We find low frequency 1/f noise that goes as L . Our values are 
in good agreement with recent experimental measurements of critical current noise in Al/AlO^/Al 
Josephson junctions. We also investigate the sensitivity of the noise on the nonuniformity of the 
tunnel barrier. 

PACS numbers: 74.25.Sv, 74.40.+k,74.50.+r 



Remarkable progress has been achieved in making 
high-quality Josephson junction qubits P, 0, d, 0, H, @, 
though sources of noise and decoherence continue to be 
problematic. Recent experiments 0, S] indicate that 
a dominant source of decoherence is two-level systems 
(TLS) in the insulating barrier of the tunnel junction as 
well as in the dielectric material used to fabricate the cir- 
cuit. These TLS fluctuators produce low frequency 1/f 



critical current noise Sj c 0, [Hi II > H, 3 ■ However, a 
simple microscopic model showing this is still missing. In 
a common scenario for Sj c fl2j , defects in the oxide tun- 
nel barrier affect nm sized conducting channels. But how 
can defects a few angstroms in width have a noticeable 
effect on a superconducting wave function with a micron 
sized coherence length that is orders of magnitude larger 
than the perturbing defect? The answer is that the tun- 
neling current is exponentially sensitive to perturbations 
of the tunnel barrier. We have confirmed this with a 
microscopic calculation of S/ c due to fluctuating TLS in 
the barrier and obtain good quantitative agreement with 
experiment. 

Previous theoretical work postulated that the qubit 
was coupled to fluctuating defects by puttin g a coupling 
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131, but no 



detects py pu ti 
term into the Hamiltonian 0, 0> Ell 2 
one has shown how this coupling arises microscopically. 
In this Letter we calculate the 1/f critical current noise 
Si c due to thermally fluctuating TLS that have electric 
dipole moments. We assume that the current I through 
the Josephson junction is given by J = I c sin 5 where I c is 
the critical current and 8 is the phase difference between 
the superconductors. In a qubit the junction is small 
(A,/ S> y/A where Xj is the Josephson penetration depth 
and A is the area of the junction) so that the phase differ- 
ence is uniform in the plane of the junction. We start by 
calculating how a dipole modifies the junction's potential 
barrier U(r). We then use a WKB formalism to com- 
pute the tunneling matrix element Tlr ~ exp VU^j 

between the left (L) and right (R) electrodes. The crit- 
ical current J r is proportional to (\Tlr\ 2 ) averaged over 
the junction [l9(. We consider elastic electron tunneling 
where different orientations of the dipole correspond to 
different values of Tlr and hence, I c . We can obtain 



Si c since each fluctuating dipole behaves as a random 
telegraph variable that has a Lorentzian noise spectrum. 
By averaging over the standard TLS distribution, and 
each dipole's orientation and position along the z— axis, 
we obtain Si c . At low frequencies we find 1/f behav- 
ior for Si c , and our values are in good ag reement with 
the corresponding experimental values 13,|13,[2l|. Our 
model predicts that the noise is very sensitive to the tun- 
nel barrier thickness L and that Sj a ~ L 5 , implying that 
the noise can be greatly reduced by decreasing L. 

In our model the TLS sits in the insulating tunnel 
barrier with an electric dipole moment p consisting of 
a pair of opposite charges separated by a distance d. 
The superconducting electrodes are located at z = 
and z = L and kept at the same potential. The an- 
gle between p and the tunneling direction [z— axis) is 9$. 
The TLS has a double well potential with a Hamilto- 
i(Aer z + AoiXr) where Aq is a tunneling 
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matrix element, A is the energy difference between the 
right and left wells, and a x ,z are the Pauli spin matrices. 
The energy eigenvalues are ±E/2 where E = \J A 2 + Aq. 
The TLS couples to the strain field. So an excited two- 
level system can decay to the ground state by emitting 
a phonon. The longitudinal relaxation rate is given by 
[22 ] Tf 1 = aEAlcoih(E/2k B T), where the prefactor a 
is a material dependent constant. The distribution of 
TLS parameters can be expressed in terms of E and T\. 
P{E,T X ) = P Q /{2T 1 ^l~T min (E)/T 1 ) [Hill where P 
is the TLS density of states. The minimum relaxation 
time T m i n (E) = T\(E = Ao) corresponds to a symmetric 
double-well potential. 

We start by calculating the junction barrier potential 
energy U when it is distorted by the presence of TLS with 
electric dipole moments 24|, |25[ . We use cylindrical co- 
ordinates (p, (f>, z) with z normal to the plane of the junc- 
tion. For a square barrier, U (p, <f>, z) = Uo — eVdi P {p, <f>, z), 
where Vdi P (p, <fi, z) is the contribution of the electric 
dipole to the barrier potential and Uq is the height of 
the unperturbed uniform square barrier. We use Green's 
functions [5^| to calculate Vdi P (r) = J G(r, r')a(r')d 3 r', 
where the charge density a(r) = Yli=i 2 1i^ 3 ( r ~ r with 
the positive charge qi = \p\/d at ri = zqZ, and the neg- 
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FIG. 1: The potential energy barrier U versus z for different 
values of p. The dipole lies along the z-axis (20 = 0.334, 
60 = 0, and d = 0.0867). We have used p = |p| = 3.7 D 
and Uo = 1 eV H§]. The inset shows the geometry and the 
parameters of the model. 



ative charge 52 = —qi at r2 = d sinc^ox + (zq + dcos6*o)z. 
The Green's function is [H G(r,r') = l/(7reL) x 
J2n=l Em=-oo sin(n7T2;/L) sn^rwr.es'/X) exp[im(0- </>')] x 
I m (nnp < /L)K m (rmp > /L), where p< = min(p, p'), p> = 
max(p, p'), 7 m and if m are the m-order modified Bessel 
functions, and e = eo£ r is the permittivity of the dielec- 
tric. The dipole potential is 



J0,Z0 j 



(1) 



where f n (p, <f>; 9 , z ) = K (mrp/L) sin(mrz /L) - 
K {nir\p—po\/L)sm{niT(zo+dcos8o)/L]. Here \p—p \ — 
(p 2 + d 2 sm 2 e - 2pdsm0 coscj)) 1 / 2 . In Fig. [fl the bar- 
rier potential U is plotted versus z = z/L for different 
values of p = p/L. The distortion of U decreases as the 
radial distance p increases, eventually becoming negligi- 
ble when p ~ 1. When the dipole flips 180°, U switches 



from (£/rj 



eVdip). 

and use the WKB approxima- 



Next we follow [2J, 
tion to calculate corresponding to the two orienta- 
tions of the dipole: 



IT, 



LR 



(p,4>)\ 2 ~exp{-2y dzyj^ ($Te^(^,z))), 

(2) 

where $ = Uq — Ek z and £k, is the energy of the electron 
incident along the z axis. We assume that eVdip <C $ 
which allows us to expand the square root in the exponent 
in powers of eVdi P /<&, and to use the WKB approximation 
which is valid for a potential that varies slowly along the 
tunneling direction. To a good approximation $ is a 
constant representing the maximum barrier height [25j ]. 
We average |7^j| 2 over the junction area A, and we find 
to lowest order in (eVdi P /<&)'- 



PF*| 2 > 



LR I 
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dppcxp|±e^/ — F(p, <£)J 
(3) 



where \T^ R \ 2 ~ exp(— 2L^j2m$/h 2 ) is the square of the 
tunneling matrix element in the absence of impurities and 
JF(p, 4>) = Jq dzVdi P (p, (j), z). T(p, (j>) will be small due to 
the oscillation in Vdi P , but it will be nonzero if there is 
asymmetry in the dipole's position along the z-axis or in 
its orientation 9q. In the exponent F(p, 4>) can still have 
a noticeable effect on the noise. Integrating over z yields 
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dp pexp[±/3 W(p, 0; 6» , zoj], 
(4) 



where the constant (3 = (p/d)/(TT 2 e) x e y /2mJ{K 2 W) and 

w(p,<p;8 ,z ) = YZ=A l - (-mfn(p,^o Q ,z )/n. 

The critical current I c is proportional to (\Tlr\ 2 ) [19|. 
Hence if = h^T^R?) I\ T lr\ 2 where h is the critical 
current in the absence of any dipoles. The critical cur- 
rent fluctuations, defined as A/ c ((9o, zq) = /+(#(), zo) — 
I~{6 0l z ), are given by 
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dp p{exp[/3 W(p, 4>; 9 , zo)} 



exp[-/3 W(p,<t>;6 ,z )]} = 



_ L 2 Ag(9 ,z ) 



A 



(5) 



where the L 2 factor comes from introducing p = p/L. 

To find the critical current noise power S/ c , we assume 
that each dipole produces a Lorentzian spectrum [27j |: 
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(6) 



where i denotes the ith dipole, P± = exp (=F£'/2fcsT) jZ 
is the Boltzmann probability of being in the upper 
(lower) state of the TLS, the partition function Z = 
2 cosh (E /2ksT), and T is temperature. We average 
over the distribution of TLS to find the low frequency 
(ujT m i n <C 1) 1// noise power: 
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where V 
tor of 1/ a/ 1 - r mm 



4/ A 

AL, Em ^> ksT, and we can neglect the fac- 
(E)/T\ because r mi „/Ti <C 1 in the 
region (E/(2k B T) > 0.1, loTi > 0.1) that dominates the 
integral. (Ag 2 ) contains the average over each dipole's 
orientation and position along the z— axis: (A5 2 ) = 

l/[2(5o M -5bJ] x /'dflosinflo/jj^ d5 A«? 2 (0 o ,5o), 
where £q lies between 5q to = d and z$ M = 1 — d to 
ensure that the dipole lies entirely in the dielectric re- 
gion. Note that the critical current noise power scales as 
L 5 /A. Although the A -1 dependence is well known ex- 
perimentally 
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111 . Il2| , it would interesting to check 



experimentally the L 5 scaling of the noise predicted by 
our model. 
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FIG. 2: Ag 2 averaged over the dipole orientation angle 6q 
versus the dipole position zq. 



To estimate (Ag 2 ), we evaluated the integrals numeri- 
cally with p m i n =0.1, which is comparable to an atomic 
radius. As Fig. Q] shows, the lower cutoff p m in is needed 
because the effect of the dipole on the tunnel barrier is no 
longer weak for p ~ 0.1. Fig. [2] shows Ag 2 averaged over 
#o versus position for various values of p ma x- The results 
overlap for p max > 2. Note that the largest contribution 
to Si c comes from dipoles near the electrodes. 

Our numerical estimates follow. We obtain (Ag 2 ) = 
1.5782 x 10~ 2 for p — 3.7 Debye (which corresponds to 
the dipole moment of an OH - impurity [2a|). d = 0.13 
nm, e r = 10, p max = 4.0, and (7 = 1 eV Q. For 
a Josephson junction with A = 1 fixa , L = 1.5 nm, 
P Q = 10 45 (Jm 3 )' 1 0, / = 1 Hz, and T = 100 mK, 
we obtain a noise power of Si c /I 2 — 4.13 x 10~ 14 Hz -1 , 
in good agreement with the recent tunnel junction resis- 
tance measurements by the Delft group (l3j who deduced 
a value of 2.0 x 10 -14 Hz -1 for the normalized critical 
current noise power using similar junction parameters. 
They also found that the noise varied linearly with tem- 
perature in agreement with Eq. jJJ . Other measurements 
in large superconducting junctions find larger S7 C than 



those of the Delft group [9|, 110, lUj- Using Eq. J7j), we 
find the critical current noise in larger junctions at 4.2 
K to be roughly 100 times lower than the experimen- 
tal values in 1CJ, 11|. The noise level difference between 
the Delft group [13j and those working with larger junc- 
tions may be due to differences in fabrication techniques. 
For example, our calculation indicates that if the effec- 
tive thickness L of the tunnel barrier in the Nb-AlO^-Nb 
junctions or Nb-NbOz-Pbln junctions d,[Tl| were 

larger by only a factor of 2 with respect to our L = 1.5 
nm, this "discrepancy" would be essentially removed. A 
T 2 dependence was observed in Ref. [llj which may im- 
ply an additional mechanism for the noise (l8j . 

The dipoles that produce critical current noise are also 
responsible for the fluctuations in the induced charges 
on the superconducting electrodes that produce charge 
noise. The dipole making 180° flips induces charge fluc- 
tuations (AQ — \2pcos9o/L\) on the electrodes. The 



time series of the charge is a two-state random telegraph 
signal with a transition rate T^ 1 resulting in a Lorentzian 
charge noise spectrum 27j ■ The charge noise results from 
averaging over the TLS density of states 17, 27|. At low 
frequencies one obtains 



^1 = -VP,k R T 
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(8) 



We estimate Sq/c 2 = 1.83 x 10 3 Hz in good agree- 
ment with the experimental values ranging from 10~ 4 
Hz" 
Eqs. 



[201 ] to approximately 4 x 10 4 Hz 1 [2l|. Using 



and ([7]) , we obtain the ratio between Sq / e 2 and 



Si a /I 2 at low frequencies: 
S Q /e 2 

sun ' 



(9) 



where B = 4/ (3(Ag 2 )). For the physical parameters 
listed above, we obtain B ~ 84.5 and (Sq/c 2 ) /(Si J I 2 ) ~ 
4.4 x 10 10 . This is consistent with the value of ~ 2 x 10 10 
deduced from experimental measurements of Sq [2l| and 
Si c [IH, though the Sq and Sj a measurements were not 
made on the same samples. 

Our calculation of the critical current noise assumed 
that in the absence of the dipoles, the tunnel barrier is 
uniform. However, local fluctuations in oxide thickness 
or barrier height make the tunnel barrier nonuniform and 
may result in "pinholes" in the barrier [29(. We can in- 
vestigate the sensitivity of the critical current noise to the 
nonuniformity of the tunnel barrier by placing the fluc- 
tuating dipole in a cylindrical island of radius p„ with 
the axis of the cylinder along z. We set the tunneling 
barrier of the island to be different from that outside the 
island: U(p) = Ui n for < p < pi n and U(p) = U ou t 
for pi n < p < pmax- We can model a pinhole by having 
Uin < U ou t- The results in Fig. [3] show that the noise 
is enhanced as the radius p~i n of the inner cylinder grows 
because the dipole distortions of the barrier potential are 
a larger fraction of £7 in than of U out - Unlike before, we do 
not expand the argument of the square root in Eq. 
When pin — > Pmax, we obtain our previous result for a 
dipole in a square 1 eV barrier. 

In the other case where t/j„ > U out , and the dipole 
is on the inner island, then the noise decreases as pi n 
increases. This is shown in Fig. [3] with EZj„ = 10 eV and 
Uout = 1 eV. The large barrier on the inner island reduces 
the amount of tunneling current through the island and 
hence limits the noise due to the fluctuating two level 
system. In the limit (pi n —> Pmin) that the inner island 
disappears, we recover our previous result for a dipole in 
a square barrier of 1 eV. This investigation shows that 
the noise power can be affected by the nonuniformity of 
the barrier, and could explain why several experimental 
groups have measured different values of Si c /I 2 [ll], 
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FIG. 3: Normalized critical current noise for a nonuniform 
tunnel barrier with two islands (Ui„ = 10 eV and U ou t = 1 
eV and vice- versa) versus the radius of the inner island pi„. 
The dashed line corresponds to noise from a dipole in a square 
barrier with a height of 1 eV. The critical current fluctuations 
in all cases are normalized by the critical current through a 
uniform square barrier of 1 eV with no dipole present. 



To conclude, we have used a microscopic model of fluc- 
tuating dipolcs in the tunnel barrier to calculate the crit- 
ical current noise in Josephson junctions. By consider- 
ing a simple tunnel barrier potential with a fluctuating 
electric dipole associated with the presence of two-level 
systems in the barrier, we estimated the critical current 
noise which was in good agreement with recent experi- 
ments. We also found that nonuniformities such as pin- 
holes in the barrier can affect the noise. It would be 
interesting to test our prediction that the critical current 
noise goes as L 5 , and to see if the ratio of the charge noise 
to critical current noise for the same junction is given by 
Eq. at low frequencies. 
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